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Abstract. We elaborate on the principle that for gapped quantum spin systems with 
local interaction “local perturbations [in the Hamiltonian] perturb locally [the ground 
state]”. This principle was established in [2], relying on the ‘spectral flow technique’ 
or ‘quasi-adiabatic continuation’ m to obtain locality estimates with sub-exponential 
decay in the distance to the spatial support of the perturbation. We use ideas of [8] 
to obtain similarly a transformation between gapped eigenvectors and their pertur¬ 
bations that is local with exponential decay. This allows to improve locality bounds 
on the effect of perturbations on the low lying states in certain gapped models with 
a unique ‘bulk ground state’ or ‘topological quantum order’. We also give some esti¬ 
mate on the exponential decay of correlations in models with impurities where some 
relevant correlations decay faster than one would naively infer from the global gap of 
the system, as one also expects in disordered systems with a localized ground state. 


1. Introduction 

This article concerns quantum spin systems with local interactions, whose dynamics 
typically satisfy a Lieb-Robinson propagation bound m- We consider paths of Hamil¬ 
tonians that are modified only locally in a confined spatial region and that maintain a 
gapped sector in the spectrum. We ask to what (spatial) extent the eigenvectors from the 
gapped sector are affected by such modifications. This setting in particular includes the 
study of ground states. The regions where the Hamiltonian is modified, may be viewed 
as impurities. 

More specifically, we investigate the principle that “local perturbations perturb lo¬ 
cally” (LPPL), which was formulated in [2]. The authors of [2] give a general and rigorous 
account of the ‘spectral flow technique’ (also named ‘quasi-adiabatic continuation’) intro¬ 
duced by Hastings in seminal work E3, see also [183, which allows to relate the spectral 
projections of an isolated patch of spectrum along the Hamiltonian path by means of a 
quasi-local unitary flow. The spectral flow is quasi-local in the sense that the involved 
unitary operator can be approximated by a truly local one up to an error that decreases 
(only) ‘sub-exponentially’ in the diameter of the support of this local operator. As an 
application, they show the principle LPPL to hold in the sense that for a unique ground 
state the effect of the perturbation’ on the expectation value of local observables decays 
(only) sub-exponentially in the distance. 

Instead of using the spectral flow technique, we relate the individual eigenvectors from 
the gapped sector (along the path) by transformations that are indeed exponentially local 
by building on ideas from [5] , based on Lieb-Robinson bounds. These transformations 
are the content of our main Theorem Id.21 in Section [3j As a drawback, these transforma¬ 
tions do not correspond to a unitary linear map between the spectral subspaces (In fact, 
they do not even correspond to a linear map, unless one gives up exponential locality). 
Moreover, they do not allow as such to consider perturbations in the whole volume, in 
contrast to the spectral flow (in that case however, it is usually hard to verify that the 
spectral gap condition remains satisfied along the path, see e.g. [26]). Here we are exclu¬ 
sively concerned with perturbations that have bounded support and hence norm-bounds 
independent of the volume. 
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By invoking the exponential clustering property for gapped ground states (see Theorem 
E2D our main result enables us to improve the LPPL principle to exponential precision 
for systems with unique ground state (also covered already in 0) or with topological 
quantum order, thus justifying our title. Furthermore, we may also dope these systems 
with impurities that increase the dimension of the low lying sector by locally closing 
the gap. This is modeled by adding on certain sites additional degrees of freedom with 
degenerate energies. In that setup, we arrive at statements on individual eigenstates 
from that sector without making further assumptions on the spectrum in that sector 
(e.g. that it contracts to a point in the thermodynamic limit, see if6|). We also give 
a result on the exponential decay of correlations in such impurity models. There, the 
decay rate is determined by the gap to the sector instead of the energy differences within 
the sector. Indeed, for disordered systems with many-body localization, one expects 
exponential decay of (most) correlations whereas the system is probably gapless. These 
applications are described in Section [4] 

Similar impurity models (restricting to perturbations of classical systems in the bulk) 
were studied before by Albanese PQ. An interesting result of this flavour was given by 
Yarotsky on ground states that are perturbations of classical states. In [25] he showed 
that the influence of arbitrary boundary conditions, which in particular may also close 
the gap, decays exponentially into the bulk, underlining the importance of a ‘bulk gap’ 
for decay of correlations besides the true spectral gap. 

Finally, in Section^ we ask the natural question whether one can construct an exponen¬ 
tially local and unitary spectral flow, or whether there is some fundamental obstruction to 
the existence to such a flow. This question is particularly important in view of the wealth 
of applications of this technique, e. g. concerning (topological) ground state phases and 
their stability Eli], the quantum Hall effect El. or disordered systems jl5], and see also 
min in. In general, there seems to be no obvious way to obtain such a exponentially 
local unitary spectral flow, but we present a simple example where it can be obtained by 
exploiting the fact that in this instance the issue reduces to a few-particle problem. 


Acknowledgements. It is a pleasure to thank Sven Bachmann for interesting discussions 
and help at the beginning of this project. Furthermore, WDR and MS are thankful to 
the DFG (German Science Foundation) and the Belgian Interuniversity Attraction Pole 
(P07/18 Dygest) for financial support. 


2. Setup 

Let r be the infinite volume, the set of vertices of an infinite graph equipped with 
the (shortest) graph distance metric d. Typically we think of the zz-dimensional lattice 
r = , but we will also need T obtained by contracting a connected subset of vertices S 

to a single vertex. We write X (g T to indicate a finite subset of the infinite volume and 
if we speak of a volume A we refer to a finite connected subset of T. For a given length 
l > 0, we call Xi := {x £ A; d(X, {x}) < /} the /-fattening of X. 

To each vertex x £ T we assign a finite dimensional Hilbert space 1~L X and to each subset 
IgT the tensor product Hx '■= The operators on T~Lx endowed with the 

operator norm constitute the algebra Ax of observables that, we say, have support in X. 
For X C X' any A £ Ax can be identified with a local operator A ® 11 acting on Rx' and 
we will not distinguish between the two in notation. Note that we do not require all 
to have the same dimension, which would be unnatural since we think of some vertices as 
corresponding to a region S C TA. 

An interaction $ on T is a family of self-adjoint local operators {$(A)}xgr, 6 

Ax- We say that it has a finite range R > 0 if <f>(Y) = 0 whenever diam(Y) > R 
(diameter of X). For each volume A the interaction $ defines a Hamiltonian (with open 
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boundary conditions) 

Ht := £ 

XCA 

and a one-parameter group of automorphisms {r t A }, t G R, on ,4 a through 
rf- A (A) = t a (A) := e iH * t Ae~ iH * t , A G A a 


( 2 . 1 ) 


( 2 . 2 ) 


which is the (Heisenberg) time evolution of the quantum spin system. The set of those 
sites within a subset IcT which interact with sites of the complement X c , 


d^X := {x € X ; 3Y <s T with x G Y and 


inYyif), rnh/0, $(F) ^ 0} 


(2.3) 


will be called the <f>-boundary of X. 


2.1. Lieb Robinson Bounds. Here we briefly state a version of the Lieb-Robinson 
bounds on propagation in quantum spin systems as proven in m- See also m by the 
same authors, where it is more explicit that the Lieb-Robinson velocity (see v below) can 
be defined independent of the single-site interaction. We fix a family of non non-increasing 
functions : [0, oo) —> (0, oo), y > 0, with F^d) := e~ /J,d F 0 (d ), such that, for y = 0, and 
hence for all /i > 0, the following hold: 

(i) uniform summability 


Fy. || := sup ^2 F r ( d ( x ’ V)) <0 ° 


xer 


y &r 


(ii) convolution property 


._ SUD V- F,{d(x lZ ))F,{d(z iy )) 

«*,„>) <0 ° 


(2.4) 


(2.5) 


To bring locality on the lattice into the game, we define a Banach space of exponentially 
decaying interactions ^ M (T), consisting of $ such that 


ii*iu ~ 


sup 

x,y£T 


E 

X3x, y 


ll*(*)ll 

F fl (d(x,y)) 


< oo 


( 2 . 6 ) 


The above expression without the contributions from single-site interaction terms, i.e., 
imposing additionally the restriction |X| > 1 in the sum, is abbreviated with ||$||( t . The 
following Theorem states the announced propagation bounds, which first appeared in [Si- 


Theorem 2.1. Let y > 0 and $ G ^ M (T). For every volume A <s T and pair of local 
observables A £ Ax and B G Ay with supports X,Y C A, the time evolution rf' A satisfies 

||[rf’ A (A),H]|| < ||A||||H||min{|5 $ A:|,|d $ r|}e-' i [ d( - Y ’ 1 ' ) - , ' |i|1 (2.7) 


for all times 1 £ K where 

211*11 ^ 

v := --- 

M 

is the so-called Lieb-Robinson velocity. 


( 2 . 8 ) 


Note that the above estimate, and in particular the Lieb-Robinson velocity v, does not 
depend on the single-site interaction terms 4>({a:}), iff, and also not on the volume A. 
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2.2. Gapped Sectors and Local Perturbations. We would like to study the effect on 
eigenfunctions of the previously defined Hamiltonians H * when adding local perturbations 
that have support only on a finite number of sites k\ , k? ,... within a subset denoted by 
K <e T. For each of these sites let Wi : [0,1] —> A{k t } be a smooth path of self-adjoint 
operators starting at Wj(0) = 0. We denote with Cw > 0 the maximal norm of the 
derivative of W := ]T7 IF,- Then 

H a (s) :=H% + W(s) (2.9) 

defines a smooth path of Hamiltonians. The restriction to (independent) perturbations 
at single sites is not severe since we are working with general infinite volume graphs and 
single-site Hilbert spaces. Generalization to perturbations involving several sites can be 
achieved by starting out with a modified infinite volume, in which these sites are merged 
to a single vertex. Note however that our estimates are sensitive to the norm of W (e.g. 
via Cw) and hence in practice to the size of K. 

Our main assumption is that H\(s) maintains a sector in the spectrum which has a 
gap of at least g > 0 to the rest of the spectrum, uniformly along the path and for all 
volumes A under consideration (e.g. for all volumes large enough to exclude non-interesting 
pathologies). Therefore the dependence on A will be mostly suppressed in notation. More 
precisely, we assume as in [2] that the spectrum of H ( s ) is a union 

c(s) = cr in (s) U CTout(s) (2-10) 

where <7i n (s) is contained in an interval of length at most L > 0 which does not intersect 
o’out(s)- We require that the distance between these two parts of the spectrum is bounded 
below by g > 0 uniformly for all s £ [0,1] (and all volumes A). 

The orthogonal projection onto the D-dimensional sector eigenspace is denoted by 
P(s). For each s £ [0,1], let A,;(s), i = 1,..., D, represent the eigenvalues in the sector of 
H(s) repeated according to their multiplicity and {^(s)}^ a set of associated orthonormal 
eigenvectors H(s)'i/ji(s) = A i(s)ipi(s). 

Given the above setting, it was shown ( 0 ) that there exists a continuous differentiable 
path of unitaries U(s) inducing the previously mentioned ‘spectral flow’ 

P(s) = U(s)P(0)U(s)* (2.11) 

and solving the evolution equation 

-id s U(s) =G{s)U(s), 17(0) = 11 (2.12) 

with the generator G(s) constructed explicitly as a self-adjoint quasi-local operator acting 
non-trivially only close to the perturbation. The accuracy with which G(s) and hence 
also 17(a) can be approximated by a local operator in Ax is shown to increase ‘sub- 
exponentially’ in d(X c ,K). Again, our main motivation here is to complement these 
findings by a result with exponential accuracy. This comes at a cost: we are able to 
transform to the perturbed eigenvectors by means of quasi-local operators which decay 
exponentially in d(X c ,K ), but this transformation will not be unitary. 

3. Main Result 

3.1. Weak Local Perturbation. In this paragraph we restrict to weak perturbations 
in the sense that we compare the eigenvectors in the sectors of H(sq + e) and H(sq) for 
so 6 [0,1 — er] and e > 0 small enough. However, the following observations are in fact the 
main step in our analysis and the result on the influence of local perturbations along the 
path H (s) will be obtained by iterating the argument. 
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Iii what follows, we assume that an interaction $ £ £$ fl ( r) is fixed, as described above, 
and that the path of Hamiltonians H(s) satisfies the main assumption stated in Section 
12.21 Therefore, constants appearing will be understood to depend only on the choice of 
p and the parameters related to the interaction $, i.e. ||Fo||, C^, and v, and also on the 
parameters involved in the main assumption: the sector dimension D , the number of 
perturbation sites \K\ and g , Cw, and L. Most importantly, constants are independent of 
the volume A. The statement ‘for e small enough’ and 0(e) is to be understood similarly, 
the maximum value can depend on all these parameters but not the volume. 

By the existence of the gap g > 0 and basic perturbation theory |9] 

\\P{s 0 + e)-P(8o)\\=0(e) (3.1) 

and the vectors {P(so + e)i/ , i(so)}i=i,...,L> remain linearly independent, if e > 0 is small 
enough. Then the perturbed eigenvectors (normalized) can be written as 


i>i(so +e) = y j c ij (so,e)P(so + s)^ j (s 0 ) 

3 


(3.2) 


for suitable coefficients Cij(so, s) £ C, i,j = 1,..., D. One rather obvious bound on these 
coefficients is 

||c,(s 0 ,e)||i := 5^|c,j(s 0 ,e)| < 2 Jd (3.3) 

3 

for e > 0 small enough. Finally we define 


c . 9+^cjn^ i 1 2 , 




(3.4) 


which will be the relevant length scale for our locality estimates. The second equality 
shows that it is related to the natural length scales in the system. 


Proposition 3.1. For every p! < l/£, there is Eq > 0 so that, for all e < Eq, sq € [0,1 — e], 
i = 1, ..., D, and l > 0, there exists a local operator P((so, e) with support in the l-fattening 
Ki and satisfying 


-Ri( s o,e) < |o-*n(s 0 + e)| 

(3.5) 

(P(s 0 + e) - R\(so, e))V’*( s o)|| = e _M ' z O(e) 

(3.6) 


Proof. In short, the proof relies on an idea from [5] to approximate spectral projection 
of gapped spectrum by a Gaussian integral of the time evolution operator, which then is 
accessible by perturbation theory and the Lieb-Robinson bounds. 

First recall that changes of the interaction at single sites as from the perturbation 
W(s) do not affect the interaction strength ||<f>||( i and in particular not the Lieb-Robinson 
velocity v. Therefore the following holds uniformly throughout the path, i.e., for every 
path parameter so■ Here we view V := H(so +e) — H(sq) as small perturbation of H(so), 
which is bounded by ||P|| < eCw- We also abbreviate Ho = H(sq), H = H(so + e) and 
the sector projections Pq = P(sq) and P = P(so+^)- For any eigenvalue A £ a = er(so+£) 
and for any given parameter a > 0 and we introduce the operator 

Vx := (ct/tt)- 1 / 2 [ d te~ at2 e lt{H ~ x) (3.7) 

JR 

By inserting the spectral decomposition 

H-\ = J2( k -\)Q k (3.8) 

where Q K is the spectral projection for H and eigenvalue re, and noting that the Fourier 
transform of Gaussian is again Gaussian, we find that 

Px = ]Te-^ K - A ) 2 Q K 

k£l<j 


(3.9) 
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Then there is a linear combination of these operators 

V:= 0 < «a < 1 (3.10) 

AGO"in 

where the sum runs over all eigenvalues in the sector a[ n = (Ti n (.so + e), that satisfies (we 
also set u"out = ^"out(^o £•)) 

V-P = J2 H e“^ (K “ A)2 g K (3.11) 

KGo-out A£(Tin 


By spectral perturbation theory ||<5 K Po|| = 0{e) for n € a ou t ; and by the gap assumption 
we therefore obtain 

\\(V - P)P 0 \\=e-^ g2 -0{e) (3.12) 

In a next step each exponential in (1X71) is replaced by its Dyson-Phillips series, which 
is norm convergent in our finite dimensional setting, 

e *t(ff-A) = (V i n [ dti ... [ n 1 d t n T°JV) ... r t 0 i (P))e l4(ff °- A) (3.13) 

n—0 do do " 7 


Here r t ° indicates the time evolution belonging to the Hamiltonian Hq and the [n = 0) 
term in the sum is taken to be the identity. If if)® is an eigenvector for an eigenvalue 
in the sector of Hq, then we can rewrite 


V$ = V a x {a/i t )~ 1/2 [ dte ~ at2 e^ H ~ x) e- lt{H °- x) e lt ^- x) ^ 
Ae<x In JR 

= V axia/n )- 1 / 2 [ dt e ~ at2 e it(x °~ x) 

A£ CTi „ Jr 

00 pt p s n — l 

J2 in dH ... / di„T t ° n (V)... T t ° (V) $ 
n—0 d0 do 

= (Rf T + R? T W 


(3.14) 


^>T < CT' 

where we also introduced a parameter T > 0 and where the operators R~ and R~ are 

defined by restricting the t-integration to \t\ > T and |t| < T respectively. The (n = 0) 

< CT’ 

term is exclusively added to i?~ . 

^>T <T 

We now show that both terms in Ri := R~ + R~ can be well approximated by a 
local operator with support in Ki if we choose 


g(g + ^m') 


dfxl 


and 


giving the three equalities 


El 

4a 


aT 2 , g(l — vT) 


4c^m'i 

■ (g + Wp ||$||> 

ygi 

g + 4C' A1 ||4>||^ 


(3.15) 


(3.16) 


which will be the dominant exponents in the next estimates. 

In fact, the effect of Rp T is so small that its locality is irrelevant. A computation gives 
the following upper bound on its norm, 


2 V a A (a/^)" 1/2 f dte~ at2 t\\V\\e tm =e~ tl ' l O(e) 

Jt>T 


(3.17) 


for any 


1 < 


ng 

g + 4C' m ||$||^ 


(3.18) 


and e small enough. Concerning the second term we can make use of the Lieb- 

Robinson bound at small times |i| < T for the time evolution of the perturbation T t °(H). 
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Let A € Ax be any operator with support X C A \Ki, i.e., at least a distance l away 
from the perturbation in K, then its commutator with Rf T is bounded by 


i[Ai?nn< 


2||Fn 


'-'LL -i '*'• 

^ n —1 

= \\A\\e-»' l O(e) 

for e small enough. Therefore (see e.g. [21 [22]) the normalized partial trace 


Rr T ' 1 := 


-Tr A \ Kl (Rf T ) 


n, ^A\K, 1 

dimft A \ ifi 

which clearly has support in Ki , satisfies 

\\Rf T ~Rf T ’ l \\=e-»' l O(e) 

By the definition of the constants a and T, equation (13.1211 becomes 

\\(V - P)P 0 \\ = e-»' 1 0(e) 

Putting it all together we find that (13.61) in the proposition holds for R\(so,e) := Rf T ’ 1 
and that (1X51) can be obtained from 


(3.19) 

(3.20) 

(3.21) 

(3.22) 


||-Ri|| < ^2 ax — l' Tin ( s o)l 

A£{Xj n 


(3.23) 


3.2. Iterating the Argument. A relation between eigenvectors ipi{ 0) and ipi( 1) from 
the beginning and end of the path H(s ) can be obtained by iterating the step (13.21) from 
the previous paragraph, 

n 

iMi) = x (n C ikik -1 ( k/n, l/n)^P{l).. .P(2/n)P{l/n)ip io (0) (3.24) 

k—1 

for n > 1 large enough (in the same sense as e small enough). In each term of the sum, 
P(m/n) then can be substituted by R l im i ((m — l)/2,1 /n) as from the Proposition. To¬ 
gether with the bound (13.31) on the coefficients Cij , the following Theorem is an immediate 
consequence of the Proposition. This is the result that motivates the title of this article. 

Theorem 3.2. For a Hamiltonian path H(s) as introduced above and for every fi' < l/£, 
there is a constant C > 0, such that, for any l > 0 and i = 1,..., D, there exist local 
operators Lr-, j = 1 with support in the l-fattening I\i, which take eigenvectors 

from the sector of H{ 0) to those of H{ 1) according to 

11^(1)(3-25) 
and which are uniformly bounded by ||Lb|| < C. 

Note again that the constant C may grow with the number of steps n needed to follow 
the path in sufficiently small steps of size e = 1 /n (see the Proposition) and therefore may 
depend on all the model parameters. 
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4. Applications 


We present two examples where the above results give a useful insight. 


4.1. Impurities in Systems with Unique Ground State. Let <f> £ -^(r) be as 
above. However, we now assume additionally that each Hamiltonian H ® has a spectral 
gap 7 > 0 above a non-degenerate ground state energy. Let if 9S be a normalized ground 
state vector, which is hence unique up to a phase factor. To each perturbation site k £ K 
we attach a finite dimensional Hilbert space I k depicting the degrees of freedom of an 
impurity. In other words we redefine the single-site Hilbert space to become ji k ®Ifc- The 
interaction <f> will be naturally viewed as family of operators on the larger space which 
act trivially on I := The energy of each impurity and its coupling to the rest of 

the system is described by a self-adjoint operator Wk on Hk ®Ik- To put ourselves in the 
context of Proposition 13.11 we assume that H = H® can be connected to 

H':= H+ W := H+ ^W k (4.1) 

k 

by a smooth path of gapped Hamiltonians satisfying the main assumption of Section 12.21 
with a dim(I) = D dimensional sector. In particular, everything holds uniformly in the 
volume A which is therefore suppressed in the notation. By basic perturbation theory, 
the main assumption is surely satisfied if the norm of W is smaller than 7 . If all I k = C 
are merely one dimensional Hilbert spaces, the overall Hilbert space and dimension of 
the ground state sector is not extended at all. Our results then effectively describe the 
influence of local perturbations on systems with unique gapped ground state, which was 
studied in detail in 0 - 

Let {ip'i}, i = 1, • ■ •, D, be an orthonormal set of eigenvectors from the sector of H'. 
Let P’ be the orthogonal sector projection for H' and P the one for the sector of H. Note 
that the range of P is the subspace of product vectors of the form ip 9S ® <j> with </> £ I. 


Proposition 4.1. For every p! < l/£, there is a constant C > 0, such that for every 
length l > 0 and local observable A with support outside of the l-fattening Ki 

\(ti,A^-(r^Ar s )\<C\d*K l/2 \\\A\\e~^ i = l,...,D (4.2) 

For every l > 0 there furthermore exists a local operator Ti with support within Ki which 
transforms the sector projection with exponential accuracy 

\\P' -T;*PT ; |[ < Ce-^' 1 (4.3) 

The norm ofTi is bounded by C. 

The Proposition can be proven by combining Theorem [372] with the exponential cluster¬ 
ing Theorem for unique ground states, which we include here for the reader’s convenience: 


Theorem 4.2 ([M]). Let p > 0 and £ & M (r) as before. Given g > 0 there exist 
constants p',c > 0 such that, if H® has a g-gapped ground state energy for a volume A 
and if Ll is a normalized ground state vector of H®, then 

|(H,AHH)| <c||A||||B||min{|a*Jf|,|a*y|}e _ '*'‘ ,( - x: ’ ir) (4.4) 


for all local observables A,B with support in X,Y, which satisfy PALI = PA*Ll = 0 with 
P being the spectral projection for the ground state subspace. The decay parameter can be 
chosen as 


g + 4C m ||<i>||( 1 


(4.5) 


This type of Theorem dates back to mmm- In case that the ground state energy 
is a non-degenerate eigenvalue of H® (as is relevant for the present proof), we can drop 
the condition on A and instead of (S3) we obtain 

| (H, ABLl) - (LI, ALl)(Ll, BLl)\ < c||A||||H|| min{|d*X|, |0 <& Y|}e _#1 ' d(x ’ i ' ) 


(4.6) 
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for all observables A, B with support in X,Y. 


The sector projection can be approximated as in (14.31) because of the simple product 
structure of vectors in the unperturbed sector, but this does not seem to be implied by 
the Theorem in general. However, recall once again that in general (14.31) holds with a 
sub-exponential bound by the spectral flow technique, see (12.111) . In physically relevant 
models \d®Kp 2 \ usually grows only polynomially in l, so that this factor in (S3 can be 
eliminated by choosing a slightly smaller decay parameter p!. 


Proof of Proposition \f.l\ First we choose an orthonormal basis {ifi} in the ground state 
subspace of the unperturbed Hamiltonian iJ of the form ipi = ip gs ® cf>i, i = 1 ,...,£>, 
with cfi G X and in Dirac’s notation we define the local operators Iij = | acting 

non-trivially on X. By Theorem 13.21 there exist, for any length l > 0, local operators Lb 
with support in Ki , so that the eigenvectors i/)' in the system coupled to impurities satisfy 


$ - 51 EVv 

j 


< (const.) e M 1 


(4.7) 


where here and in the following ‘(const.)’ stands for a constant that may change from 
line to line, but may only depend on the fixed model parameters (not the volume). The 
greatest of these possibly different constants qualifies as the constant C in the Proposition. 
p! is of course taken from the Theorem. Therefore the perturbed sector projection can be 
approximated according to 


pi _ \ ' tI t pj tI* 
M / j 1 J ip- L pi ± J iq- LJ iq 


< (const.) e M 1 


(4.8) 


and we arrive at (14.31) with 




Concerning the first part of the Proposition note that 

M'i) - - ( const -) e -//i/2 


(4.9) 


(4.10) 


again with the aid of Theorem 13.21 By the exponential clustering Theorem (Theorem 03 

| (Vv> AI mp Ll /p 2 *L l !q “ (V’m. A i>m) (V’m) 1 L 'lq 1 q rJ>m) \ , 4 u . 

< (const.) \d®Ki/ 2 \\\A\\ e _M 1 


for any m = 1,... ,D and a p' > 0 as specified in the Proposition, which finishes the 
proof. 

□ 


4.1.1. Exponential Decay of Correlations. Just above we stated and used the well-known 
fact that gapped unique ground states arising from local interactions exhibit exponential 
decay of spatial (truncated) correlations. The same was shown to hold for states defined 
through spectral projections on gapped sectors of eigenvalues which approach the ground 
state energy in the thermodynamic limit nans],- For our model with a single impurity we 
show a (modified) exponential clustering property for every state defined by an eigenvector 
from the low-lying sector. In particular it holds for any ground state (in finite volume) 
independent of the gap to the energy of the next highest excitation, which may be small 
but non-vanishing for all volumes. This is not surprising since we just showed that each 
sector state is identical to a gapped unique ground state away from the impurity and since, 
intuitively, one would expect the correlation length to be determined by bulk properties, 
e.g. a kind of ‘bulk gap’, of the system. 
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For two subsets X, Y C A, which will be the supports of observables A and B, and 
depending on the perturbation set K, we define an effective distance dx{X,Y) as the 
smallest number l > 0 such that X; U Y; U Ki has a connected component containing both 
X and Y. This distance hence tends to be shorter if an impurity is located in between 
the supports X and Y. Our next result shows that correlations decay exponentially in 
this effective distance for a single impurity K = {/c} and we assume that \d®Ki\ increases 
at most polynomially in l. 

Proposition 4.3. If \K\ = 1 and for every p! < l/£, there is a constant C > 0, so that, 
for every pair of observables A, B with support in X , Y 

< C'||A||||B||e-^ x ’ r ) (4.12) 

for i = 1,... ,D. 


Proof. If d(X, Y) is smaller than or equal to the distance of X and Y to the perturbation 
set K, then exponential clustering for w' is an immediate consequence of Proposition lj.il 
knowing that lo 9S has this property, see Theorem 14.21 More generally, we use again that 
i’i ~ E j L ij J ji ip i as in (EE2D and that Lb commutes with either A and/or B , let’s say B , 
for l = dx(X,Y)/ 2. In the same way that (14.101) was obtained in the previous proof we 
arrive at 


(4>i, ABfti) ~ 


T T AT l T Ftih 

m i 1 mp J - / ip-™" L 'iq 1 qm 1 -* ti 


J < (const.)e ^ 1 (4-13) 


p,q 


Since the distance d(X U Ki,Y) between the support of each I mp L l * p AL\ q I qm and of B is 
at least dx{X,Y) and since we assumed that \d®Ki\ grows slower than exponentially in 
l. the exponential clustering Theorem can be used to finish the proof. 


□ 


For several impurities \K\ > 1, the different regions around them may be correlated 
for arbitrary distances for states from (entangled) eigenvectors already due to possible 
degeneracies in the spectrum. This can of course already be seen in the uncoupled case 
with W = 0. For the average of sector states w(-) = tr(P •) however, truncated corre¬ 
lations still decay exponentially in the distance between observables as a consequence of 
the exponential clustering Theorem applied to the unique ground state in the bulk. If the 
impurities can be coupled to the system one after another along a gapped path, one might 
also expect exponential clustering to hold for u/(-) = tr(P' •) as above in Proposition 
m (i.e. with respect to distance dx), but we are not able to show it. 

Also to clarify this statement, we briefly point out that in this case ‘sub-exponential’ 
decay of correlations can be easily obtained with the help of the spectral flow technique, 
indicating the deficiency in Theorem 13.21 of not being implemented by a unitary transfor¬ 
mation. The set of impurity sites can be split into two parts K = K x U K y of sites close 
to X and close to Y in the sense that, for l = dx(X, Y)/2, their fattened sets satisfy 

d(K x U X, Kj U Y) > / (4.14) 

Furthermore we assume that both the uncoupled and perturbed Hamiltonian H and H’ 
are connected to 

H x :=H+ Yl W k (4-15) 

kGK x 

by a gapped path in the same way as before, but with derivatives in and A^y re¬ 

spectively. This requirement formalizes the assumption that the impurities can be coupled 
to the system one after another without closing the gap. The spectral flow technique [2] 
allows to construct two unitary operators U x and U y supported on K x and K} which 
transform between the gapped sectors of H/H x and H x /H\ 

U?P{U?)*kP x and UlP x {U^)*^P' 


(4.16) 
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up to errors whose norm decays sub-exponentially in l indicated by the ss symbol. Here, 
P x stands for the sector projection of H x . Then indeed 

u'(AB) « u((U x )*{U X )*ABU X U x ) 

= u{(u x yAu x (u x yBu x ) 

« u{{u x yAu x )u({uYyBu x ) ( 4 . 17 ) 

= u({u x y{u?yAu?u x )u J ((u x y(jj x yBuru x ) 

^uj'(A)uj , (B) 

To get to the third line, the exponential clustering property was used for the state w. 

4.2. Impurities in Systems with Topological Quantum Order. Another interesting 
class of models for which Theorem l3.21 translates into exponentially sharp locality estimates 
are those with topological quantum order (defined below). To account for the topological 
aspects of such systems, each volume A here is a finite graph imbedded on a possibly 
non-trivial surface rather than a finite subgraph of a common infinite volume T. For the 
sake of concreteness, we restrict to two dimensional square lattices A = Z^xZiof length 
L with periodic boundary conditions (imbeddings on a torus) and to translation invariant 
finite-range interactions. This restriction includes one of the most prominent examples 
with topological order, Kitaev’s toric code, see m ■ To relate the models at different 
L , we fix a p and we assume that the model parameters || Fo||, C M , and v converge to 
finite non-zero values as L —> oo. This puts us in a setting where Theorem 13.21 applies 
(A-uniformly). 

We are concerned with impurity models which are set up for each such A in the same 
way as in the previous section and the same notation and assumptions will be used. The 
only change we make is that we dismiss the requirement of a unique ground state for H 
and instead we assume, for each volume A, the ground states of H possess topological 
quantum order (TQO) defined as in HSUS], see also p?B1 : 

For any local observable A whose support X is contained within a square of side length 
L* there is z G C, so that PAP = zP. 

This implies that different ground states cannot be distinguished by measurements local¬ 
ized on length scales up to L*. Typically we have in mind that L* increases as a function 
of the system’s size L , e. g. that L* > L a for some a > 0. 

Proposition 4.4. For every p! < l/£, there is a constant C > 0, such that the following 
holds for every length l > 0: If A is a local observable with support X outside the l- 
fattening Ki of the impurity set K such that the union X U _K)/ 2 is contained in a square 
of side length L*, then 

K^A^-zS^KCWA] |e-^ (4.18) 

for i,j = l,...,D and where z £ C is given by PAP = zP. 

The Proposition shows that away from the impurities in the system (or from perturba¬ 
tions if all Ik = C) the measurements of local (up to L*) observables remain unchanged 
with exponential accuracy. Moreover, in this region TQO persist ‘with error’ that is ex¬ 
ponentially decreasing. This terminology was coined in [F;. 

Proof. We essentially follow the proof for Proposition 14.11 with adaptations due to the 
possible, say /-fold, degeneracy of the ground state energy for the bulk system (without 
the impurities). Let {ifi} be an orthonormal basis in the ground state subspace for H of 
the form 

A = Vf* ® </>i K . * = (ib, Ik) 

*6 = 1,...,/, i K = l,...,dim(X) 


(4.19) 


12 


W. DE ROECK AND M. SCHUTZ 


Note that D = f ■ dim(I). By using Theorem 13.21 and the local transition operators hx-.m 
as before we find 


y 'X ^(vh.mx)'’ ^mxVK^iy ^jq ^qKm K ^(qb,mx)) 
p,q 

< (const.) e _M 1 / 2 


(4.20) 


for any rriK = 1,..., dim(I). By the definition of the local transformation operators and 
assumptions of the Proposition, the support of 

AT rl/2*r!/2. 

■^^rriKPK ip jq qnrrix 

lies within a square of length L* and therefore those terms in the sum with pb ^ qb vanish 
as a consequence of the TQO condition. We can assume that the support X of A is disjoint 
from K otherwise there is nothing to be proven. On the square lattice |<9 <l> A'// 2 | only grows 
polynomially in l. Again by the TQO condition 

i^(pb,m K ) ’ ~~ "fipbqb ( 4 - 21 ) 

and in particular 

P(A - z)Tjj [pbimK) = P(A* - z)i> {pbimK) = 0 (4.22) 

for all pb and ttik, so that the proof of the Proposition follows by the exponential clustering 
Theorem 14.21 

□ 


5. Simple Example of an Exponentially Local Spectral Flow 

Here we present a particular impurity model in the same setup as above for which it 
is not difficult to construct a unitary exponentially local spectral flow. We are dealing 
with a p-dimensional lattice of S = 1/2 spins, i.e. T~L X = C , and restrict to cubic volumes 
A = Z£ of length L with periodic boundary conditions. The bulk interaction is that of an 
‘xy-model’ with Hamiltonian 

H = - Y, (SlSl + SlSD+YW^ + ^^/Z + Sl) (5.1) 

d(x,y)=l x 

where S'*, i = 1,2,3, are the Pauli matrices at site x. For u we take a positive function 
bounded below by u(x ) > 7 , for 7 > 0 , which then is also the spectral gap above the zero 
ground state energy. The ground state is the product of ‘spin-down’ vectors (eigenvectors 
of S x and eigenvalue —1/2) and frustration-free. Note that (1/2 + S x ) is the orthogonal 
projection onto the ‘spin-up’ state and that in the first sum we can replace S x Sy + S x S p 
by S+S- in terms of the spin raising/lowering operators S± = S x ± iS x . In fact, this 
model is unitarily equivalent to a system of hard core bosons where sites occupied by a 
boson correspond to those being in spin-up state (sometimes called Matsubara-Matsueda 
correspondence BSD- Locality of operators is strictly preserved in this correspondence. 
Using standard second quantization notation, see e.g. |7J, 

H^H-= P hc dT(-A + u) P hc (5.2) 

where Phc is the orthogonal projection on the subspace of bosonic Fock space T(/ 2 (A)) 
with at most one Boson per site (hard core condition). A denotes the discrete Laplacian. 
The particle number is conserved in this many-boson system just as our xy-model con¬ 
serves the spin, i.e. [H, = 0. 

We now add an impurity at a single site, K = {k}, which itself consists of N spins 
with Hilbert space I = (C 2 )® Ar . This is equivalent to adding N sites I = {i \,... An} 
on the Boson model side of the correspondence. The coupling to the system is described 
by a smooth path W{s) of spin-conserving operators on Hk with H( 0) = H, so that 
H(s) = H + W{s) maintains a 2 Ar -dimensional y-gapped sector. Recall that we denoted 
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with Cw the maximal norm of the derivative of W. Initially, any state with Bosons located 
only at the Impurity’s sites is a ground state. By the assumption of spin/particle number 
conservation along the path, the low-lying sector is fully described by eigenfunctions of a 
few-particle system (up to TV). These decay exponentially, making the model so tractable. 


Proposition 5.1. Let P(s) be the spectral sector projections for the low-lying sector along 
the path of Hamiltonians H(s) in the xy-model with impurity as defined above. There are 
C, p 1 > 0 depending only on g and TV, such that the following holds: for every length 
l £ [0,oo], there is a smooth path of self-adjoint operators Gi(s) and of unitary operators 
Ui(s) with support in the l-fattening Ki of the impurity site, which satisfy 

||P(s) - f^;(s)P(0)f/*(s)|| < Ce~»' 1 (5.3) 


and solve 

-id a U l {s) = G l (s)U l (s), 17,(0) = H (5.4) 


We only give a brief sketch of the proof: We show that Kato’s ‘transformation function’, 
see II. § 4 in [5], which is a particular choice of a spectral flow, is exponentially quasi¬ 
local for our model. It is defined as the unique solution U(s) of (15.411 generated by the 
commutator 


G(s) :=i[P(s),d 8 P{s)] (5.5) 

and satisfies (15.31) with l = oo. For n > 0, let H( n \s) and p( n \s) be the restrictions of 
H ( s ) and P(s) to the invariant n-spin (particle) subspace. If n > TV, note that P(s ) = 0, 
which implies G^ n \s) = 0 and U^ n \s) = 11 for the restrictions of G(s ) and U(s). By the 
holomorphic functional calculus, the sector projection and its derivative can be expressed 
in terms of the resolvent as 


d s p( n) (s) = -J- ( d z{H^(a)-z)~ X d a W^{a){H^(s)-z)~ l 
27r * Jc{s) 


(5.6) 


for a family of contours C(s) at distance of at least g/2 to the spectrum and enclosing the 
low-lying spectrum. Switching to the particle description of our model, (s) is unitarily 
equivalent to to an n-particle discrete Schrodinger operator on (the symmetric hard-core 
subspace of) / 2 ((AUJ)“). For such operators, a Combes-Thomas type estimate [B] shows 
that the matrix elements of the resolvent as in the above integrals decay exponentially 
away from the diagonal in the canonical position basis. The length scale of this decay 
may increase with the dimension n and as the gap closes (for z £ C(s)). Since d s Pt n \s) is 
equivalent to an operator with support in KUI, the matrix elements of G^ ( s ) then decay 
exponentially away from k. Truncating G^ (s) outside of K\ only gives a difference in 
norm that decays exponentially in l. Therefore, we can define Gi(s) of the Proposition as 
the direct sum (from 0 to TV) of these truncations and Ui(s) as the corresponding unique 
solution of differential equation (15.41) . 
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